On the flexoviscous behavior in viscoelastic materials 
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A phenomenological model on the relationship between viscoelasticity and flexoelectricity is pro- 
posed to address the electromechanical coupling under rheological processes in viscoelastic materials. 
Our theoretical studies show that some observed ferroelectric phenomena in certain viscoelastic ma- 
terials, like polymer films, may not be intrinsic material behavior. 

PACS numbers; 



There is currently a great deal of interest in electrome- 
chanical coupling phenomena in polymer films due to 
their potential technological applications. Although the 
piezoelectric effects (both the direct and the converse 
ones) and the corresponding pyroelectric effects in poly- 
mers had been observed before 1970s, the existence of the 
ferroelectricity in polymers was not recognized until the 
end of 1970s (Ref. [ij and references therein). However, 
there has been a concern about whether these phenom- 
ena are intrinsic material behavior or improper proper- 
ties due to the inherent rheological or viscoelastic proper- 
ties of polymers. Recent studies on the flexoelectricity in 
polymers have indirectly shown that the electromechani- 
cal coupling in viscoelastic materials may involve certain 
non-trivial collective and cooperative material behavior, 
which may not be explained by using conventional piezo- 
electric and flexoelectric models proposed for solid crys- 
talline materials [H-Q- Those experimental observations 
stimulate us to consider the electromechanical coupling 
and the ferroelectric behavior in viscoelastic materials. 
To address our problems clearly, we start our discussion 
from a brief survey of piezoelectric and flexoelectric ef- 
fects. 

For non-centrosymmetric solid crystalline materials, 
their piezoelectric effects can be written as ^ 

Dk = (1) 

and 

S^i ^ '^ijk^k, (2) 

where Dk and Ek represent the induced electric displace- 
ment and the applied electric field, respectively; Tij and 
Sij represent the applied stress and the induced strain, 
respectively. Eq. ([T]) and Eq. ^ are usually called the 
direct piezoelectric effect and the converse piezoelectric 
effect, respectively; d^-^. is defined as the piezoelectric 
coefHcient for the direct piezoelectric effect and dj^j, the 
piezoelectric coefficient for the converse piezoelectric ef- 
fect. For Hookean materials, we have df^f. = d~-j. as 
required by equilibrium thermodynamics. 



The inhomogeneous electromechanical energy coupling 
between the applied strain gradient and the induced elec- 
tric polarization in solid crystalline materials, i.e., the di- 
rect type II flexoelectric effect, was proposed by Kogan 
in 1964 0, which is given below. 

^™-M-^feV,5.„ (3) 

where is the symbol representing the gradient with 
respect to the direction k; Mi/mfc defined as the flexo- 
electric coefficient for the direct type II flexoelectric ef- 
fect. For Hookean materials, we can easily write, via the 
Maxwell relations, the converse type II flexoelectric effect 
as 0,i: 

= l^ijmk'^kEjn, (4) 

whereuf/~i, is defined as the flexoelectric coefficient for 
the converse type II flexoelectric effect. Similarly, for 
Hookean materials, we have fJ-lj^n^ ^ t^ljmk- 
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FIG. 1: (a) Uniform strain field; (b) gradient field of electric 
field; (c) gradient field of strain; (d) uniform electric field. 
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In 2007, Fu and Cross considered a situation where 
the Landau-Ginzburg-Devonshire (LGD) phenomenolog- 
ical model may not work properly in a gradient field of the 
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order parameter in ferroelectric materials and proposed 
that there are two types of flexoelectric effects, the type I 
and the type II, existing in all solid crystalline dielectrics 
Q . The phenomenological model of the type II flexoelec- 
tric effects is indeed the Kogan's phenomenological model 
given in Eqs. jS]) and (j4]), whereas the type I flexoelectric 
effects are novel inhomogeneous electromechanical cou- 
pling phenomena predicted by Fu and Cross. The con- 
stitutive equations of this phenomenological model can 
be written as [1]: 

VkD^m - ^^l;„,,s.J (5) 

and 

"^kTij^ til-^.^f^Era, (6) 

where and ^[^^4, are defined as the flexoelectric 

coefficients for the direct and the converse type I flexo- 
electric effects, respectively. Clearly, M^jt^fe = ^J^lJmk ^^"^ 
Hookean materials. 
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FIG. 2: A piece of material undergoing shearing driven by a 
force F; (a) F = 0; (b) F / 0. 



The physical pictures of the flexoelectric effects can 
be diagrammatically explained in a two-dimensional case 
shown in Fig. 1; mechanical energy and electrical en- 
ergy are uniformly distributed in Figs. 1(a) and 1(d); 
electrical energy and mechanical energy are distributed 
in certain vortex patterns to form the electric displace- 
ment gradient shown in Fig. 1(b) and the stress gradient 
shown in Fig. 1(c). The energy coupling from Fig. 1(a) 
to Fig. 1(b) and from Fig. 1(d) to Fig. 1(c) represent 
the direct and the converse type I flexoelectric effects, 
respectively; the ones from Fig. 1(c) to Fig. 1(d) and 
from Fig. 1(b) to Fig. 1(a) represent the direct and 
the converse type II flexoelectric effects, respectively. In 
both cases shown in Figs. 1(b) and 1(c), the magnitude 
of the vertical electric displacement values and the verti- 
cal stress values reaches the maximum at spots 1 and 3 



but these values are oriented in opposite directions, and 
reaches zero at spots 2 and 4; we can see the same sit- 
uations at spots 2 and 4 and then at spots 1 and 3 in 
the horizontal values. From the viewpoint of equilibrium 
thermodynamics, these gradient fields are energetically 
more favorable since the vector values with opposite di- 
rections would cancel each other so that the resultant 
Gibbs free energy will decrease and approach a minimum 
value. Therefore, for a material undergoing such inhomo- 
geneous electromechanical energy coupling, its molecules 
or atoms will be forced to move or rotate by a flexo- 
electric torque, which is generated during the coupling 
process by the inhomogeneously coupled energy, to form 
the aforementioned gradient fields. Unfortunately, com- 
pared with crystal binding energies, such a torque is too 
small to change crystalline structures of solid dielectrics 
so that the aforementioned vortex structures or gradi- 
ent fields cannot be formed under normal conditions in 
solid dielectrics. However, for viscous materials like liq- 
uid crystals that can deform even under a very small 
stress, flexoelectric effects, especially the type I flexo- 
electric effects, could significantly change their molecular 
alignments and alter their properties [9[; for viscoclastic 
materials, flexoelectric effects may also have an influence 
on their properties. To make a reasonable and mean- 
ingful connection between flexoelectricity and viscosity 
in viscoelastic materials, we re- visit the concept of shear 
viscosity in purely viscous materials. 




(a) (b) 



FIG. 3: (a) a perfect crystal (we use a zig-zag pattern to 
diagrammatically represent the crystalline structure); (b) a 
partially crystallized polymer with both the crystalline struc- 
ture and the amorphous structure. 



We consider a viscous material undergoing a rectilinear 
shearing flow, in which only the x-componcnt of the rate 
of shear deformation is nonzero and is a linear function 
of y, driven by a force as diagrammatically shown in Fig. 
2; its shear viscosity can be simply written as [loj : 

Tx = r]dux/dy, (7) 

where is the shear stress in the x-axis direction; 
represents the rate of shear deformation in the x-axis 
direction; 77 is defined as the shear viscosity. For the vis- 
cous material, there is no relationship between the shear 
deformation and the driving force. 
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(c) 



(d) 



FIG. 4: Diagrammatic sketch of the formation of the tilted 
structures in polymers; (a) a polymer with a vertical molec- 
ular network; (b) the tilted structures induced by both flexo- 
electricity and flexoviscosity in the polymer given in (a); (c) 
a polymer with a horizontal molecular network; (d) the tilted 
structures induced by both fiexoelectricity and flexoviscosity 
in the polymer given in (c). 



For a Hookean material also diagrammatically shown 
in Fig. 2, Eq. ([7]) does not hold. The relationship be- 
tween the shear deformation and the driving force can be 
defined as t = G9 {6 <^ 1), here 9 = a/h describes the 
shear deformation; t = F/A {A is the top surface area 
of the material); and G is the shear modulus. There- 
fore, when F is removed, the shear deformation can be 
recovered completely. 

For a viscoelastic material diagrammatically shown in 



Fig. 2, when F is removed, the shear deformation can be 
recovered to a large extent shortly after the occurrence of 
the deformation; however, after longer times the recov- 
erable deformation reduces. This is because the material 
possesses both viscous and elastic properties simultane- 
ously. So we can re-write the Eq. (O as: 



7]d{Sxy)/dt, 



(8) 



where Sxy is the strain representing the gradient of the 
shear deformation in the y-axis direction. A universal 
shear viscosity relationship in viscoelastic materials is 
given below 11 



i]ijkid{Ski)/dt, 



(9) 



where i]ijki represents the usual viscosity tensor. 

Substituting the strain tensor term in Eq. ^ into 
Eq. ([S]), we obtain the following relationship between 
piezoelectricity and viscosity in non-ccntrosymmetric vis- 
coelastic materials, 

Tij = r]ijkid^i^^d{Em)/dt = -fijmd{Em) / dt, (10) 

where "fijm is defined as the piezoviscous coefficient. This 
equation demonstrates a linear relationship between the 
applied time- varying electric field and the induced shear 
stress. The phenomenon is defined as the piezoviscous 
effect or piezoviscosity. 

By taking the partial derivative of both sides of Eq. © 
with respect to m, we obtain a modified equation shown 
below, 



V 



r]ijkid(yjriSki)/dt. 



(11) 



Substituting the strain gradient term embodied in Eq. 
^ into Eq. pT|) . we have the following relationship be- 
tween the type I fiexoelectricity and viscosity in viscoelas- 
tic materials, 



J 



VfcT, 



kTij ~ rjij 



.^d{Era)/dt 
I 



Pijmk 

d{Em)/dt 



(12) 



where /9f,,„jj. is defined as the type I flexoviscous coeffi- 
represents the compliance tensor. The equa- 



cient; s 



IJUV 



tion describes a linear relationship between the applied 
time-varying electric field and the induced inhomoge- 
neous shear stress. For simplicity, we define this phe- 



nomenon as the type I fiexoviscous effect. 

Similarly, by substituting the strain tensor term em- 
bodied in Eq. (|4]) into Eq. ([9]), we have the following 
relationship between the type II fiexoelectricity and vis- 
cosity in viscoelastic materials. 



J 



= VijuvS^Juv^J-"„,kdi'^kE,n)/dt = pj^^^^^diV k E„,) / dt , 



(13) 



where pfj,^k defined as the type II flexoviscous coef- ficient. Eq. ([T^ demonstrates a linear relationship be- 
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twccn the applied time-varying inhomogeneous electric 
field and the induced shear stress. The phenomenon is 
defined as the type II flexoviscous effect. Both type I 
and II flexoviscous effects, for simplicity, can be called 
flexoviscosity. 

It is interesting to note that c?^;^ appears in Eq. 
(Uni), and ^J-lf^k /if/m/c appear in Eq. (HI]) and Eq. 
(HSl), respectively; but and cannot 

be found in those equations. Of course, d^^^, A*[j>nfcj 
and fJ-lj^^. should have certain relationships with the vis- 
cosity tensor, but these relationship would be different 
from Eq. dTO]), Eq. dE]), and Eq. which means that 

coelastic materials. 

Perhaps, it is more interesting to use the concepts 
of flcxoelectricity and flexoviscosity to analyze dielectric 
and mechanical behavior of polymer materials since they 
are not only viscoelastic materials but also possess the 
following two characteristics: (a) electric charges can be 
easily trapped and stored in polymers, especially some 
highly insulating polymers, for extremely long periods of 
time; therefore, they are often regarded as quasi-electret 
materials jl^. Thus for polymer materials, Em given in 
Eq. g]), Eq. dSl), Eq. (Uni), Eq. and Eq. 1^ should 

be substituted by E^* + E^, here E!^* is the applied 
electric field and E^^ is the electric field generated by the 
trapped charges. This means that, even without E^*, 
the fluctuation of E^ due to the thermally excited hop- 
ping of electrons from one isolated state to another in 
polymers could gradually alter their long chain molecu- 
lar structures via flcxoelectricity (the short-term effect) 
and flexoviscosity (the long-term effect); (b) the molecu- 
lar structures of polymers arc usually very complicated. 
A perfect crystal and a semicrystalline polymer are dia- 
grammatically compared in Fig. 3. Being different from 
a perfect crystal, a semicrystalline polymer has inher- 
ent heterogeneous structures like a ^^spontaneous com- 
posite^^ formed with crystalline parts in a glassy matrix 
p^ . The mutual interactions between crystalline phases 
and amorphous phases and among different amorphous 
phases in polymers could dominate their material prop- 
erties. For example, piezoelectric responses have been 
detected in non-piezoelectric polymer fllms 0, 0] . There- 
fore, the electromechanical coupling in polymers should 
be certain collective and cooperative material behavior, 
which is different from that in solid piezoelectric materi- 
als. 

From Fig. 3(b), we can see that neither electrical en- 
ergy nor mechanical energy can be uniformly distributed 
in polymer materials. So it is convenient to use the con- 
verse type I flcxoelectric effect deflncd by Eq. ([6]) and 
the type I flexoviscous effect deflncd in Eq. ([T2|) to an- 
alyze the formation of flexoelectric domains or gradient 
fleld related structures in polymers. Here, we only give a 
qualitative analysis on the formation of these structures. 
We assume that there is a single polymer layer with a 
simple vertical molecular network shown in Fig. 4(a); the 



network would be altered via flcxoelectricity and flexo- 
viscosity to form a tilted molecular structure shown in 
Fig. 4(b); similarly, we can see another layer with a hori- 
zontal molecular network shown in Fig. 4(c) and the cor- 
responding layer with a tilted molecular structure shown 
in Fig. 4(d). Imagining that a polymer film consists of 
many single layers shown in Fig. 4, we can see that some 
tilted and layered structures would be formed and pre- 
served (due to viscoelasticity) in the film; the ensemble of 
these structures looks like a helix-like structure in three- 
dimensional space. Before discussing the influence of 
these structures on material properties of the fllm, we re- 
visit Meyer's pioneer work on ferroelectric liquid crystals. 
In 1970s, Meyer proposed that the liquid crystal with the 
chiral smectic C phase or the smectic C* phase should be 
ferroelectric [iJl- Meyer's creative argument is based on 
the point group symmetry analysis; he argued that chi- 
ral molecules lack the inversion symmetry and their tilted 
and layered phases would be non-centrosymmetric, which 
renders the existence of ferroelectricity in the liquid crys- 
tal. In polymers, molecules might not be chiral and the 
aforementioned tilted and layered structures might not 
be non-centrosymmetric; however, those tilted and lay- 
ered structures could create a non-uniform spatial distri- 
bution of the induced charges (induced by, for example, 
the electromechanical coupling) and the trapped charges 
and generate "a quasi-permanent electrical potential ori- 
ented along a specific direction!^ (which is different from 
the electrical potentials randomly induced by the trapped 
charges) in polymers, which could exhibit a ^^spontaneous 
electric polarization^^ that can be reversed by the appli- 
cation of an external electric field. The reason that such 
an electrical potential can be preserved over relatively 
long periods of time and docs not dissipate is due to the 
viscoelastic behavior in polymers. Therefore, we can ar- 
gue that some observed P-E hysteresis loops in polymer 
films might be improper ferroelectricity. Experimental 
studies are now underway to verify the arguments ad- 
dressed here. 

Last but not least, the process of the formation of the 
tilted and layered structures diagrammatically demon- 
strated by Fig. 4 should be accompanied by material 
property variations, especially mechanical and dielectric 
properties, in polymer films, which can be seen clearly 
in Fig. 4. This means that the phenomenological model, 
fiexoclectricity and flexoviscosity, could provide us with a 
broader perspective for understanding the physical aging 
phenomena of polymers. 

In summary, we propose here a phenomenological 
model on the relationship between viscoelasticity and 
the electromechanical coupling in viscoelastic materi- 
als, i.e., the flexoviscous effect and the piezoviscous ef- 
fect. According to the model, the thermodynamic equiv- 
alence between the direct and the converse flexoelectric 
and piezoelectric effects is broken in viscoelastic materi- 
als. The extended physical picture of the model is that 
the long chain molecular structures of polymers could 
be altered via fiexoclectricity and flexoviscosity to form 
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certain tilted and layered structures, which could in- 
duce ferroelectricity-like behavior in polymers. There- 
fore, some observed ferroelectric phenomena in polymer 
films might be improper properties. 
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